Stoichiometry, structure, and transport in the quasi-one-dimensional metal, 

Li . 9 Mo 6 Oi 7 
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A correlation between lattice parameters, oxygen composition, and the thermoelectric and Hall 
coefficients is presented for single-crystal Lio.9MoeOi7, a quasi-one-dimensional (Q1D) metallic com- 
pound. The possibility that this compound is a compensated metal is discussed in light of a sub- 
stantial variability observed in the literature for these transport coefficients. 



Introduction. Lio.gMogOi7 known as "lithium pur- 
ple bronze" (LiPB), is a low-temperature superconductor 
(T c «2K) first synthesized and studied in the 1980s.— ~— 
It has attracted interest more recently for its quasi-one 
dimensionality and Luttinger-liquid candidacy^ - — It is 
distinguished among quasi-one-dimensional (Q1D) com- 
pounds by the absence of a conventional density-wave 
transitio n 1 ^ 11 (either charge or spin) throughout a broad 
temperature range, T > T c . An upturn in its resistiv- 
ity below Tm ~ 30 K may be associated with localiza- 
tion, dimensional crossover or the development of un- 
conventional (e.g., electronically-driven) charge density- 
wave orderi 10 ' 11 

Values reported for the chain- axis (6-axis) electrical re- 
sistivity of LiPB 2 ' 10 ' 13 ^ vary by more than an order of 
magnitude, from 0.4 mf2cm to more than 10 mficm at 
300 K, and appear to be responsible for similar differ- 
ences in reported values for the anisotropy ratio, p c j Pb 
(p c is the resistivity along the next most conducting di- 
rection) . Some of the differences may be associated with 
inadequate shorting^ for current along the chains, of 
the voltage drops in directions transverse to the current 
flow, a delicate matter in highly anisotropic conductors. 
Increases in pi, for crystals^ annealed in air at 200°C 
suggest that a variable oxygen stoichiometry also con- 
tributes to the reported differences. Recent thermoelec- 
tric and Hall measurements 1 - support such a view, in- 
dicating small differences in structure and stoichiometry 
that correlate with transport. This is an important issue 
given the renewed attention being paid to LiPB and the 
sensitivity of conduction in such low-dimensional materi- 
als to defects or impurities that may alter the low-energy 
electronic structure. A variable stoichiometry that allows 
for control of charge-carrier density could prove useful in 
the study of LiPB which is proposed as a model, bulk 
system exhibiting Luttinger-liquid physics. 

The principal finding reported here is a correla- 
tion between oxygen content and the c lattice param- 
eter for a number of as-grown LiPB crystals. See- 
beck (thcrmopower) and Hall coefficient measurements 
on LiPB crystals from our own work 19 and those of 
otherS ) 15 ' 16 ' 20 i 21 are then presented together. Measured 
under open circuit conditions, the thcrmopower is not 
prone to the same challenges of transverse shorting that 
may complicate comparison of resistivity measurements. 



Hall measurements are also less susceptible to such errors 
given proper averaging with reversed current and mag- 
netic field directions. The low-T chain-axis thermopower 
varies considerably among the crystals measured?^ and 
correlates with their c-axis lattice parameters. Results 
for the Hall coefficient show considerable variability but 
the role of stoichiometry is less definitive as there are not 
enough data available on individual crystals for which 
structure, stoichiometry and transport coefficients have 
all been measured. The possibility that LiPB is a com- 
pensated metal, motivated by the behavior of the ther- 
moelectric coefficients^ is examined in light of the Hall 
data. 

Crystal Structure and Stoichiometry. Single-crystal 
growth of Lio.gMo60i7 using a temperature-gradient flux 
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FIG. 1. (color online) Examples of high-angle extrapolations 
of XRD reflections used to determine the a and c lattice pa- 
rameters for Lio.gMo60i7 single crystals. Solid lines are least- 
squares fits. Inset: rocking curves of the (500) reflections for 
two crystals. 



2 



90.64 



<D 90.62 
T3 




29 m 



9.47 9.48 9.49 9.50 

c lattice constant (A) 

FIG. 2. (color online) Lower panel: Oxygen content vs. c lat- 
tice constant for Lio.9MoeOi7 crystals. Right (left) ordinate 
scales correspond to EDS (EDS+2.5 Wt.%) oxygen contents 
(see text). Error bars represent the standard deviation of 
EDS oxygen contents from multiple locations and of high- 
angle c-parameter extrapolations from multiple planes. Ar- 
rows labeled Al, A2, and C indicate the c lattice parameters 
for specimens discussed in the text. Upper panel: f3 vs. c for 
all crystals studied. 



method is described in detail elsewhere^yLL The crystals 
grow as thin platelets, with a ~ 0.05 — 0.20 mm and 
0.5 — 2 mm for b and c. A Philips X'Pert x-ray diffrac- 
tometer (Cu K a , A = 1.54056 A) was employed in de- 
termining the crystallographic structure. Lattice param- 
eters for the monoclinic unit cell22 were determined by 
from high-angle extrapolation 2 ^ of reflections from vari- 
ous lattice planes (Fig. [IJ , and angles /3 (between the a 
and c axes) from differences in d-spacings for the (hOh) 
and (hOh) reflections. Some crystals exhibit resolution- 
limited rocking curve widths, FWHM < 0.05° (without 
a monochromator) , while others (typically larger crys- 
tals) exhibit a modest mosaicity, FWHM ~ 0.1 - 0.2° 
(inset, Fig. [T]). For the latter, rocking curves exhibit 
a convolution of two or more resolution-limited peaks, 
indicative of separate domains misaligned by small an- 
gles. A total of ten crystals were studied and their a and 
b lattice parameters were found to be the same within 
uncertainties, 12.752(2) A and 5.520(2) A, respectively. 
These agree with those determined from recent neutron 
scattering studies 24 on a powder specimen ground from 
thousands of similarly prepared crystals from our group. 

The c lattice parameters were found to vary among 
the as-prepared crystals, correlating with oxygen con- 
tent for seven of the crystals that were examined with 
x-ray energy-dispersive spectroscopy (EDS) on a sys- 
tem equipped with a polymer window suitable for light- 



element detection. Such variations, observed in the 
present study for crystals from the same growth batch 
and from different batches, is common in flux-grown crys- 
tals. The results from three or more EDS scans (with area 
1 /Ira x 1 /im) from different locations on the surface of 
each crystal were averaged to determine the oxygen con- 
tent. Li, not detectable in these scans, was not included 
in the compositional normalization [chemical analysis of 
the neutron scattering powder specimen 2 ^ indicated a Li 
content 0.924(9)]. To ensure that the crystal surfaces 
(probed to a depth of several £tm by EDS) were repre- 
sentative of the bulk, pieces of two crystals were polished 
to half their original thickness (corresponding to a depth 
of 50-100 /im) and examined by EDS for comparison. 
The average oxygen contents determined from measure- 
ments at several locations on the interior surfaces of these 
crystals differed from those of the original surfaces by no 
more than variations between different locations on each 
surface. 

Fig.[2]shows the correlation between c and oxygen con- 
tent for these crystals. Standardless oxygen determina- 
tions from EDS typically show small systematic discrep- 
ancies due to inadequate accounting for absorption. As a 
measure of this discrepancy we use the neutron scattering 
results 24 which indicate full oxygenation (17 O/f.u., cor- 
responding to 32.09 Wt% ignoring Li) for c = 9.4909(2) 
A. Comparing our data for the same c values suggests 
that the EDS oxygen results are « 2.5 Wt% too low. 
The right (left) ordinates of Fig. [^correspond to the EDS 
(corrected: EDS+2.5 Wt%) oxygen contents. 

The maximum variation in oxygen content for the crys- 
tals examined, ~ 2 Wt%, corresponds to ~ 1 oxygen 
atom per f.u. The corrected EDS oxygen scaling in Fig. [5] 
implies that some crystals have excess oxygen, presum- 
ably accommodated as intcrstitials. On the other hand, 
the fact that four of the crystals have the same mea- 
sured oxygen content within uncertainties but differing c 
parameters suggests that other factors may influence c, 
e.g. variations in the Li concentration^ or inhomogene- 
ity of the oxygen distribution (on a scale smaller than 
the 1 /xm). Though prior work 2 suggests a very narrow 
range of Li nonstoichiometry (0.87 — 0.93 per f.u.), small 
Li concentration variations between crystals may exist 
and might correlate with those observed here for oxy- 
gen. The c parameter is thus as a surrogate indicator for 
stoichiometric variation in oxygen and possibly Li. We 
note that a crystal with oxygen excess (labeled "Al" in 
Fig. [5]) was remeasured after 18 months stored in a dessi- 
cator and found to have larger c (labeled "A2" in Fig. [2J , 
suggesting a tendency of crystals with oxygen excess to 
slowly lose oxygen or for oxygen to redistribute, via dif- 
fusion, over long periods of time. Small increases in the 
angle j3 with increasing c parameter are shown in Fig. [3] 

Interstitial oxygen in other oxides (e.g. cuprates) may 
affect the electronic structure through local charge dop- 
ing, lattice distortion, or magnetic moment formation. 
The quasi-one-dimensional, conducting Mo-0 double 
chains of LiPB may be especially susceptible to perturba- 
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FIG. 3. (color online) (a) Fermi surface from Ref. |g, (b) 
schematic of the conducting Mo-O-Mo zig-zag chains along 
the crystallographic b axis, and (c) the arrangement of dou- 
ble chains within the ac plane. 



tions from such effects. The LiPB band structure^— i£ 
is composed of two nearly-degenerate, Q1D bands cross- 
ing the Fermi energy (with kp ~ it /2b), yielding two 
pairs of slightly warped Fermi surface (FS) sheets in the 
b*-c* planes with very little dispersion along a* [Fig. [3] 
(a)]. The small splitting predicted for the bands near 
Ep from density functional theory^ (~ 30 — 40 meV) 
has not been observed in photoemissioni^ The splitting 
is associated with an opposite warping of the FS sheets, 
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attributable within tight-binding parameterization: 
opposite signs (and comparable magnitudes) for the hop- 
ping integrals between constituent chains of a ladder 
and between ladders (t±2) hi neighboring unit cells 
[Fig. H] (b), (c)]. These transverse hopping parameters 
may be sensitive to variations in the c parameter. 

Thermopower. The 6-axis thermopower of LiPB, 
S b (T) (Fig. 0}, was first reported in Ref. [U More 
recently^ Sb(T) was reported for several crystals hav- 
ing different c parameters (letter labels in Fig. @| . Data 
for crystal D have been extended in the present work 
to 385 K. The inset of Fig. [4] shows Sb as a function 
of c. While the room-temperature thermopower is sim- 
ilar for all crystals, the dramatic low-T peak in Sb near 
T = 50 K correlates with the c parameter. The very large 
value (200 /xV/K) found for specimen Al with the small- 
est c parameter and highest oxygen content (Fig. [2]) is 
unusual for a metal. The largest Nernst coefficients (v), 
with maxima near T = 20 K, were founcU^ for sped- 




100 200 300 

T(K) 



400 



FIG. 4. (color online) S b (T) for crystals from Ref. Ql (cir- 
cles with letter labels) and Ref . [13 (diamonds). Inset: St at 
T =50 K and 300 K vs. c lattice parameter for crystals from 
Ref. |H. 



mens C and D having the smallest values for Sb- This 
favors a two-carrier (ambipolar) conduction scenario in 
which very large hole and electron partial thermopowers, 
S h and S e (S e < 0), nearly cancel due to compensation: 
Sb oc S h + S e , v oc fi(S h - S e ) (n is the carrier mobility). 
A departure from the compensation condition in speci- 
mens with larger c parameter could explain the appear- 
ance of their large low-T values for S b - The mechanism 
underlying such very large partial thermopowers at low 
T remains to be established, but phonon-drag appears to 
be a plausible candidate^ 

Though all specimens have positive 300 K thermopow- 
ers, their T dependencies approach a linear-in-T form 
with negative slope that results in negative values for 
Sb above room temperature (specimen D, Fig. U). This 
implies a negative (electron-like) carrier diffusion ther- 
mopower along the chains that is independent of the low- 
T behavior. 

Hall Coefficient. The Hall coefficients of crystals A2 
and C having c = 9.482(4) and 9.491(4), respectively (ar- 
rows in Fig. [2]) , were reported elsewhere Their oxygen 
contents were not measured; their Hall data are shown 
in Fig. [S] Also shown are data from three other sources 
in the literature i 15 i 16 ' 21 All of the data show very simi- 
lar T dependencies, but magnitudes vary by a factor of 
10. Note that Rh > throughout the T range for all 
measurements. 

At present it is not possible to draw conclusions about 
a dependence of Rh on c parameter with complete data 
sets only for specimens A2 and C (c was not reported 
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FIG. 5. (color online) Hall coefficient for crystals A2 and C 
(Ref. 03) an d others from the literature. Data for Ref. Il6l 
were computed by interpolation from their reported a xy /B, 
pb and p c data as, Rh = {<J xy / B)p b p c . 



along with the other data reproduced in Fig. [3]). In 
a more recent workii studying crystals from the same 
source as that of Ref. [lU, lattice parameters a = 12.73 A 
and c = 9.51 A where reported, both differing signifi- 
cantly from those reported here. Thus in addition to a 
stoichiometric variability observed for the present crys- 
tals, there may well be other differences between spec- 
imens prepared by different groups that have yet to be 
identified. 

The Hall data might be reconciled with the compensa- 
tion scenario motivated by the thermoelectric data. Em- 
pirically, that Rh is always positive (a sign opposite to 
that of the diffusion thermopower) and varies substan- 
tially in magnitude among those crystals measured tends 
to support such a view. To examine this issue in more 
detail, first note that for Q1D systems with open Fermi 
surfaces the relation between the low-field Hall coeffi- 
cient and carrier density differs from the free-electron 
expression, and for constant mean-free-path is general- 
ized to^ Rh = v/( n \ e \)i where r] measures the non- 
linearity of the dispersion along the conducting chains, 
r] = — {hkp / vp)d 2 e / dk"l | ep , and vf = de/dkb\ £F . 

From the perspective of correlation effects, LiPB is 
one-quarter filled (approximately half an electron per 
band crossing Ep) with a tight-binding dispersion, 7,8 
e(fcb) = — t cos(fc;,&/2). The Hall coefficient correspond- 
ing to k F = 7r/(2&) of LiPBi would then be^2- R H = 
-(7r/4)/(ra|e|) ~ -1.7 x 1(T 9 m 3 /C assuming 1.9e - /unit 
cell as dictated by the chemistry and bonding. Though 
this value is in the middle of the experimental magni- 
tudes near room temperature, the sign is opposite to 
those observed. The Hall coefficient for coupled Lut- 



tinger chains differs from the noninteracting case only 
by small power-law-in T corrections^ Theoretical treat- 
ments of the thermopower for single Luttinger chains 
with impurity scattering 29 indicate a linear temperature 
dependence for the carrier diffusion contribution similar 
to the noninteracting case. 

A positive sign for Rh can arise from negative curva- 
ture in the dispersion (r] > 0) along one or both of the 
FS sheets (e.g. the density functional bands 5 -^ suggest 
this is the case for the band that crosses at higher mo- 
mentum along the P — K direction), or from a strong 
momentum dependence for the mean-free path^i due 
to electron-electro n 25 ' 31 or electron-phonon interactions, 
though such scattering anisotropy tends to be weaker at 
high T. Two-carrier fitsi 9 - to the Rh data for crystals 
A2 and C yield hole and electron partial Hall coefficients 
that are an order of magnitude or more larger than Rh, 
implying carrier densities substantially smaller than ex- 
pected. A reduced mobile carrier density is suggested by 
observations of localized charge in optical studies^ as 
might be anticipated for carriers confined to chain frag- 
ments isolated by disorder. 

An alternate, intriguing possibility is that strong corre- 
lations modify the FS from that predicted within density 
functional theory [Fig. [3] (a)] . The Mott gap (A) in LiPB 
appears to be quite small 8 - such that t± > A. In such 
systems having strong on-site and longer-range Coulomb 
repulsion, calculations revea l 32 ' 33 a regime intermediate 
between ID insulator and 2D metal where the suppres- 
sion of the Mott gap via transverse hopping leads to the 
formation of a FS broken into narrow electron and hole 
pockets. This physics was recently invoked as a possible 
explanation for Fermi arcs found in photocmission stud- 
ies of Y-124 due to the CuO chains M> Such a picture is 
appealing for LiPB as it offers a possible explanation for 
two-carrier physics apparent in the thermoelectric coeffi- 
cients, positive values for Rh (provided there is a higher 
weighting of hole- like pockets), and carrier densities in- 
ferred from a two-carrier analysis of Rh that are much 
smaller than expected. Whether this scenario can be con- 
sistent with photocmission observations of the FS awaits 
a more complete and higher-resolution mapping through- 
out the Brillouin zone^ 

Summary. Variations in the c-lattice parameter of 
Lio.gMogOi7 crystals, attributable to stoichiometric vari- 
ations in oxygen and possibly Li, correlate with the low- 
T thermopower. The implication is that the low-energy 
electronic structure is sensitive to small variations in sto- 
ichiometry in as-grown crystals. Thus the c parameter 
is revealed by this work to be an important metric with 
which to compare specimens. The data suggest forma- 
tion of interstitial oxygen during growth as a possible 
source of oxygen variability. The available Hall data are 
less conclusive as to trends with c, but the observed vari- 
ations in magnitude may be consistent with a two-carrier 
picture motivated by thermoelectric data. 
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